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Abstract 

We propose an experiment, which would ahow to pinpoint the role of spin- 
orbit coupling in the metal-nonmetal transition observed in a number of two- 
dimensional systems at low densities. Namely, we demonstrate that in a 
parallel magnetic field the interplay between the spin-orbit coupling and the 
Zeeman splitting leads to a characteristic anisotropy of resistivity with respect 
to the direction of the in-plane magnetic field. Though our analytic calculation 
is done in the deeply insulating regime, the anisotropy is expected to persist 
far beyond that regime. 
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In a recent papeiS an interesting experimental observation was reported. It was demon- 
strated that the period of beats of the Shubnikov-de Haas oscillations in a two-dimensional 
hole system is strongly correlated with the zero-magnetic-field temperature dependence of 
the resistivity. The beats of the Shubnikov-de Haas oscillations have their origin in the 
splitting of the spin subbands in a zero magnetic field,iii. The authors of Ref. |I| were able 
to tune the zero- field splitting by changing the gate voltage. They observed that, while in 
the absence of the subband splitting, the zero-magnetic-field resistivity was temperature- 
independent below T = 0.7K , a pronounced rise (by 5 percent) in resistivity with tempera- 
ture emerged in the interval 0.2K < T < 0.7K at the maximal subband splitting, indicating 
a metallic-like behavior. This close correlation suggests that it is a mechanism causing the 
spin subband splitting that plays an important role in the crossover from the metallic-like to 
the insulating-like temperature dependence of resistivity with decreasing carrier density (the 
metal-nonmetal transition). This transition by now has been experimentally observed in a 
number of different two-dimensional electronS^ and holei^ systems. By challenging the 
commonly accepted concepts, it has attracted a lot of theoretical interest and attempts to 
identify the underlying mechanism. Possible relevance of zero-field splitting to the transition 
was first conjectured in Ref. |TB|. The evidence presented in Ref. |I| about the importance of 
the subband splitting for metallic-like behavior of resistivity is further supported by the very 



recent data reported in Ref. 14 



Another important feature of the metal-nonmetal transition, which might also provide a 
clue for the understanding of its origin, is that the metallic phase is destroyed by a relatively 
weak parallel mag netic field0'§-§. At the same time, no quenching of the metallic phase 
in a parallel magnetic field was observed in SiGe hole ga^Hll, in which the strain, caused by 
the lattice mismatch, splits the light and heavy holes. 

As far as the theory is concerned, the role of the parallel magnetic field was previously 
accounted for exclusively through the Zeeman energy, which either alters the exchange in- 
teractions (and, thus, electron-ion binding energyBil) or suppresses the liquid phase,@ 



or 



affects the transmittancy of the point contact between the phase-coherent region: 
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It is appealing to combine the observations§'0 of the subband sphtting in zero field and 
the resultsi'iHil in a parallel magnetic field within a single picture. The spin-orbit (SO) 
coupling appears to be a promising candidate for such a unifying mechanism. Indeed, on 
one hand, it is known to lead to spin subband splitting. On the other hand, a parallel 
magnetic field, though not affecting the orbital in-plane motion, destroys the SO coupling 
and, thus, suppresses the intersubband transitions. Possible importance of these transitions 
was emphasized in Ref. Their suppression with increasing magnetic field is caused by the 
fact that the corresponding subband wave functions become orthogonal for all wave vectors. 

At the present moment there is no consensus in the literature about the role of the 
SO coupling. Several authorS^^l have explored the role of the SO coupling as a possible 
source of the metallic-like behavior, by considering noninteracting two-dimensional system 
and including the SO terms into the calculation of the weak-localization corrections. At 
the same time, the majority of theoretical workS^'frU, stimulated by the experimental 
observation of the transition, disregarded the SO coupling. 

To pinpoint the role of the SO coupling in the metal-nonmetal transition, it seems im- 
portant to find a qualitative effect which exists only in the presence of the SO coupling. 
Such an effect is proposed in the present paper. We show that an interplay between the SO 
coupling and the Zeeman splitting gives rise to a characteristic anisotropy of resistivity with 
respect to the direction of the parallel magnetic field. Obviously, the Zeeman splitting alone 
cannot induce any anisotropy. To demonstrate the effect, we consider the deeply insulating 
regime, where the physical picture of transport is transparent. 

We choose the simplest form for the spin-orbit Hamiltonian@'@ 

ifso = ak- (cr X z). (1) 

Here a is the SO coupling constant, k is the wave vector, z is the unit vector normal to the 
2D plane, cr = (ai, (72, era) are the Pauli matrices. In the presence of the parallel megnetic 
field, the single particle Hamiltonian can be written as 
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(2) 



where m is the effective mass, g and are the g-factor and the Bohr magneton respectively; 

= 2gfiBB is the Zeeman sphtting; 0b and 0k are, correspondingly, the azimuthal angles 
of magnetic field B (Fig. 1) and the wave vector k. The energy spectrum of the Hamiltonian 
Eq. (M) is given by 
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E±(k) 



± -JaI + 4a2A;2 + AakAz sin( 
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Note that the spectrum is anisotropic only if both Az and a are nonzero. 

The standard procedure for the calculation of the hopping conductance is the following^. 
Denote with P12 the hopping probability between the localized states 1 and 2. The logarithm 
of P12 represents the sum of two terms 
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where the first term originates from the activation; is the activation energy0 and T is 
the temperature. The second term in Eq. (|^) describes the overlap of the wave functions 
of the localized states centered at points Ri and R2, so that R = Ri — R2. In Eq. @ we 
use the fact that within the prefactor the overlap integral coincides with the Green function 
G(R). For the matrix Hamiltonian Eq. (|^), the Green function is also a matrix 



G(R) 



(27r)2 ^ - ^(k) ■ 
By projecting onto the eigen-space of Hamiltonian Eq. ( 
presented as 
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the above expression can be 
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where the projection operators Pt(k) are defined as 



n(k) 



1 0(k) 
0*(k) 1 



p_(k) = l-p+(k). 



(7) 



where 0*(k) is the complex conjugate of 0(k), which is defined as 

^exp(-i0B) -m/i;exp(-i0k) 

£4k)-i5-(k) ■ 

When the distance R is much larger than the localization radius, ao, the integral over 0k 
is determined by a narrow interval |(/)k — 0r| ~ {kR)~^^^ <S 1. This allows to replace 0k by 



in the square brackets and perform the angular integration. Then we obtain 



271 r°° dkVk , 



(9) 



The next step of the integration is also standard. Namely, for large R, the G(R) is de- 
termined by the poles of the integrand. However, in the case under the consideration, the 
equation E±(k) — E leads to a fourth-order algebraic equation. To simplify the calculations 
we will restrict ourselves to the strongly localized regime ^ ma ^/h^. In this case the 
poles can be found by the successive approximations. In the zero-order approximation, we 
get the standard result k — iko, where ko is defined as 



V 2m E 

ko = = ^ ■ (10) 

In the first order approximation, we have k — iko + k\ where ki is given by 



ma 



A;i = ±i-^^A2-l + 2iAisin((/)B-0R), (11) 
where the dimensionless Zeeman splitting Ai is defined as 

Ai = i^z/'^ak^- (12) 
Within this approximation, the long-distance asymptotics of the Green function is 

G{R) oc e-^/<'t>-^^'t'^\ (13) 
where the decay length is given by 

1 - -^Re^jAj - 1 + 2iAi sin(0B - ^r) j • (14) 



In the last equation it is assumed that the real part, i?e (...), has a positive sign. Our main 
observation is that the decay length and, concommitantly, the probability of hopping are 
anisotropic, when the parallel megnetic field and the SO coupling are present simultaneously. 
By evaluating the real part in Eq. (0) we obtain 



a(0B,0R)-' = A;o(l - ^I^V^i - 1 + \/l + At-2A?cos2(0B-0R)j. (15) 

To characterize the anisotropy quantitatively, we introduce the perpendicular decay 
length a_L = a(0B — 0r = ±f ) and the parallel decay length a\\ = a((/)B = 0r). Then 
a quantitative measure of the anisotropy can be defined as 

a_L — ail ma 



where the function f{x) is given by 



■/(Ai), (16) 



/(x) =a;-(x2-l)V20(x-l), (17) 

where 6{x) is the step-function. In the strongly localized regime {ako <^ \E\) the anisotropy 
is weak. The magnetic field dependence of the anisotropy is shown in Fig. 2. It can be 
seen that the maximal anisotropy corresponds to Ai = 1 and it vanishes both in strong 
and weak magnetic fields. The theory of hopping transport in the systems with anisotropic 
localization radius is presented in Ref. The principal outcome of this theory is that 
the anisotropy of the localization radius (and, consequently, the exponential anisotropy of 
the hopping probability (^)) does not lead to the exponential anisotropy of the hopping 
resistance. In fact, the exponent of the resistance is the same as for the isotropic hopping 
with localization radius ^ya\\a±. However, the anisotropy in the Green function manifests 
itself in the prefactor of the hopping resistance^ 

^^ = C^^-C^/(Aa (18) 

where C ~ 1 is the numerical factor, determined by the perturbation theory in the method 
of invariants for random bond percolation problemill. The exact value of the constant C 
depends on the regime of hopping (nearest-neighbor or variable-range hopping) . 



The microscopic origin of the SO Hamiltonaian Eq. (HD is the asymmery of the confine- 
ment potential. In III-V semiconductor quantum wells there exists another mechanism of 
the SO coupling, which originates from the absence of the inversion symmetry in the bulk 
(the Dresselhaus mechanism@). Within this mechanism Hso = Pio'xkx — (^yky) (for [001] 
growth direction). Then the calculation similar to the above leads to the following result for 
the anisotropic decay length 



a{<pB, = ko(l- -^^]JAI - 1 + ^1 + A| + 2Ai cos(0b + 0r) j, (19) 

where A2 is related to the Zeeman splitting as 

A2 = Az/2f3ko. (20) 



By using Eq. (|T9D we get for anisotropy 

a± — a|| m(3 



^/(A2), (21) 



where the function / is determined by Eq. (|T] 

In conclusion, we have demonstrated that, due to the SO coupling, the rotation of an 
in-plane magnetic field with respect to the direction of current should lead to a characteristic 
angular variation of resistivity with a period tt. The anisotropy is maximal for intermediate 
magnetic fields and vanishes in the weak and the strong-field limits. In the strongly localized 
regime, considered in the present paper, the magnitude of anisotropy is small. However, as 
seen from Eqs. (|16D and ([T8|), the magnitude of anisotropy should increase as the Fermi 
level moves up with increasing carrier concentration (since decreases). So the resistivity 
is expected to remain anisotropic, perhaps with a modified angular dependence, far beyond 
the deeply insulating regime. For high enough concentrations the SO coupling (subband 
splitting) is negligible, so that the anisotropy should be also weak. If the intersubband 
scattering governs the metal-nonmetal transition, then the resistivity anisotropy should reach 
maximum around the critical density. 

Finally, let us discuss two possible complications for the experimental observation of the 
anisotropy in resistivity. Both of them stem from the fact that a realistic two-dimensional 



system has a finite thickness. Firstly, with finite thickness, even a small deviation of the 
magnetic field direction from the in-plane position would cause a certain anisotropy even 
without SO coupling. However, in this case, the anisotropy would only increase with in- 
creasing magnetic field, while the SO-induced anisotropy should vanish in the strong-field 
limit. The second effect of the finite thickness is that it causes the anisotropy of the Dres- 
selhaus term with respect to the crystalline axes. As it is shown in Ref. the interplay 
of anisotropic Dresselhaus and isotropic Bychkov-Rashba terms results in the crystalline 
anisotropy of the resistivity in the weak-localization regime. This effect should be distin- 
guished from the anisotropy with respect to the direction of current predicted in the present 
paper. 
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FIGURES 

FIG. 1. Azimuthal position of the in-plane magnetic field, ^b, and of the wave vector of 

electron, (pi^, are shown schematically. 



FIG. 2. Magnetic field dependence of the anisotropy of the decay length. 
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